Study of gauge (in) dependence of monopole dynamics 
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We investigated the gauge (in)dependence of the confinement mechanism due to monopole condensation in 
SU(2) lattice QCD by various abelian projections. We found (1) the string tension can be reproduced by monopoles 
alone also in Polyakov gauge and (2) the behaviors of the Polyakov loop at the critical temperature seem to be 
explained by the uniformity breaking of the monopole currents in every gauge. 
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1. Introduction 

Many people believe that the quark confine- 
ment mechanism in QCD can be understood as 
dual Meissner effect due to abelian monopole 
condensation. And this picture is supported by 
recent Monte-Carlo simulations of abelian pro- 
jected QCD in maximally abelian gauge (MA 
gauge). There are infinite ways of abelian projec- 
tion extracting an abelian component. If abelian 
monopole condensation is really the confinement 
mechanism in QCD, it must be independent of 
the way of abelian projection, i.e., gauge inde- 
pendent. The purpose of this report is to study 
monopole physics in gauges other than the MA 
gauge such as Polyakov gauge. 

The method of our studies are the followings. 
(1) We perform usual Monte-Carlo simulations 
of QCD. (2) Choosing a gauge, we extract an 
abelian theory by abelian projection [1, 2]. (3) We 
measure observables written in terms of abelian 
link variable. When the features of QCD can be 
reproduced by such an abelian variable, we call 
this property abelian dominance. (4) The oper- 
ators of the abelian observables ( Wilson loop, 
Polyakov loop, (ipip), . . . ) can be decomposed into 
contributions from monopole and from photon. 
We measure these two parts separately. When the 
features of QCD can be explained by the contri- 
bution from monopole alone, it is called monopole 
dominance. 

Many results in the MA gauge and a few results 
in the other gauges about abelian dominance and 
monopole dominance have been reported. 
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2. Study in Polyakov gauge 

We consider finite temperature QCD. The 
abelian dominance of the quantities derived from 
the Polyakov loops is trivial in the case of the 
Polyakov gauge. 

We define a matrix P(s) as follows 
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The Polyakov gauge is defined by diagonaliz- 
ing P(s). And an abelian link variable u(s,/i) 
is extracted as U(s,/i) = c(s, p)u(s, ji), where 
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u(s,/i) is diagonal. Then all U{s,4) are diag- 
onal and abelian. U(s,4) = m(s,4). It means 
that the full Polyakov loop is equal to the abelian 
Polyakov loop. The Fadeev-Popov determinant 
of the Polyakov gauge is given in a simple form 

[11]- 

The abelian Polyakov loop operator: 

N 4 

P = exp(i M s + (*' " l)4)0 4 (s + (i ~ 1)4) (2) 

8 = 1 

U(s,4) = «(s,4) = dmg{e is ^ s \ e - iSi ^} 

can be decomposed into a photon part and a 
Dirac string part using following identities: 

e 4 (.s) = -J2d( s - s')[d' v e vi {s') + s 4 ( W*'))]> (3) 

s' 

where 9^ = d^By — 8^6^. D(s — s') is the lattice 
Coulomb propagator. The current conservation 
94/4(5) = leads us to [8] 

P = Pi (photon part) • P2 (Dirac string part), (4) 



JV 4 



P 1 = eX p{-i^ J 4 (s + (i- 1)4) 

8 = 1 

j2D( s + (i-i)l- s ')d'J^( s ')}, 

s' 

P 2 = exp{-27ri J i( s + (* - I) 4 ) 

8 = 1 

yD(s+(i-l)i- 8 ')dln vi (s')}, 



(5) 



(6) 



where 

0^(s) = 0^(s) + Zirn^is), (-7T < 6^(s) < ir). 

Here *nfj, v {s) = \ifj,vpan pa (s) is the Dirac string. 
The monopole currents k^s) are the boundary of 
the Dirac string: = d v *n vfi [12]. 

The string tension is derived from the correla- 
tion of Polyakov loops. The simulations are per- 
formed on a 24 3 x 4 lattice. The measurements 
are done every 50 sweeps after 2000 sweeps for 
thermalization. We used 100 configurations for 
the Polyakov loop and 500 or 1000 configurations 
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Figure 1. Monopole (Dirac string) and photon 
contributions to Polyakov loops in the Polyakov 
gauge 



for Polyakov loop correlation. The Polyakov loop 
is plotted in Fig. 1, and a typical Polyakov loop 
correlation (/? = 2.27) is shown in Fig. 2. The 
string tension was obtained from Polyakov loop 
correlation by least square fit assuming that the 
static quark anti-quark potential is given by linear 
+ Coulomb + constant terms. The [3 dependence 
of the string tension is shown in Fig. 3. 

We can see in Fig. 1 that the abelian flux is 
squeezed by the Dirac string also in the Polyakov 
gauge. 

In the Polyakov gauge, we found that the 
abelian Polyakov loop correlations give the same 
string tension as the full one, (abelian domi- 
nance) and that the Polyakov loop correlations 
from monopoles reproduce much the same string 
tension near the scaling region, (monopole dom- 
inance) 

3. Monopole dynamics at f3 c 

Monopole dominance of Polyakov loops was 
shown in various gauges [8]. It suggests that the 
monopole dynamics which makes these behaviors 
of Polyakov loop at f3 c is gauge independent. 

Considering Ji{s) = 1 along the Polyakov line, 
the monopole contribution to Polyakov loop can 
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Figure 2. Polyakov loop correlations in the 
Polyakov gauge 
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Figure 3. [3 dependence of the string tension in 
the Polyakov gauge 



be rewritten by 

(P) = (exp(-2«^ D 3 (s - r) l -e V A P od' v V(*0)>- 

s' 

(7) 

Here D%{s) is the 3-dimensional Coulomb propa- 
gator and *fi^ v is the Dirac string projected in the 
3-dimensional space, 

*n^ v {s) = ^e^pa n po{s)- (8) 
This reduces to 

(P) = (exp(2«^)>, (9) 

47T 

where £l(s) is the solid angle made by monopole 
loops projected in the 3-dimensional space from 
the Polyakov loop (s). 

Notice that this equation shows that the value 
of Polyakov loop depends only on monopole cur- 
rents and does not depend on a form of Dirac 
sheet. 

The behaviors of the monopole loops were stud- 
ied in SU(2) QCD in the MA gauge [13]. In the 
confinement phase, there are one long connected 
monopole loop which distributes uniformly and 
some short loops. The length of the long loop 
becomes shorter as [3 becomes larger. In decon- 
finement phase, a long loop disappears. 



The properties of Polyakov loop can be un- 
derstood by these monopole dynamics qualita- 
tively. In the confinement phase (T < T c ), long 
monopole loops distribute uniformly. The SI can 
be random from to Ait. The value of the 
Polyakov loop at each point is random. The av- 
erage of the Polyakov loop is zero. On the other 
hand, in the deconfinement phase (T > T c ), there 
is a space where no monopole exists. In such a 
space, £1 takes small value and the local Polyakov 
loops are nearly one. The average of the Polyakov 
loop becomes non-zero. 

This picture seems to be correct in the MA 
gauge. In other gauges, the monopoles are so 
dense that we can not see the breaking of unifor- 
mity clearly. 

Hence to check if this picture is correct in other 
gauges, we investigated the histogram of the solid 
angle. Monte-Carlo simulations were performed 
on 24 3 x 4 lattices at j3 = 2.2, 2.3 and 2.4, j3 = 2.3 
is critical [3. We took 50 configurations 

In Fig. 4, the shape of the histogram changes 
at f3 c in every gauge. It is almost flat in the con- 
finement phase, whereas it has a peek near zero 
in the deconfinement phase. 

Notice that, in the Polyakov gauge, there are 
several methods of diagonalizing Yl Si ^4- If we di- 
agonalize it with a condition [Q > 0, the his- 
togram becomes antisymmetric. So we used the 



3 



. MA gauge A 

i i 


— P=2.2 




Polyakov gauge 


----- p=2-3 






I i 


— p=2.4 






I ' 

i 

i 




10000 




i 

j 
i 




5000 




/ 






■ ■ (5=2.2 S 








0-2.3 








— (5=2-4 






n 







£1/2 



- F 12 gauge 














— P=2.2 




|S 2 3 " 




— P=2.4 



-It 

1 0000 ■ 





n/2 



- No gauge fixing 



p=2.2 
(3=2.3 
R=2.4 





n/2 





n/2 



Figure 4. Histograms of £1/2 on 24 3 x 4 lattices 
with 50 configurations. 



method of diagonalizing it without such a condi- 
tion. 



4. Summary and Discussion 

We found the following results. (1) The abelian 
dominance of the quantities derived from the 
Polyakov loops is exact in the Polyakov gauge. 
(2) The string tension and the characteristic 
behaviors of Polyakov loops are reproduced by 
monopoles also in the Polyakov gauge. (3) The 
properties of Polyakov loops have a close relation 
to the uniformity breaking of the monopole loops 
at T c , and this picture seems to be gauge invari- 
ant. 

In Fig. 3 some difference between the values of 
the sting tension from monopoles and full one are 
seen at the strong coupling region in the Polyakov 
gauge. The study of finite size scaling in the T = 
system is important and in progress. We must 
measure the string tension using Wilson loops on 
larger lattice. 

In the Polyakov gauge, there are several ways of 
diagonalizing Yl Si ^4> anc ^ the monopole currents 
seem to depend on this procedure. We must study 
if these dependence affect physical quantities. 



In the study of the relation of the monopole 
and the string tension in the MA gauge [7], it was 
found that not all monopole currents contribute 
to the string tension. Only a long monopole loop 
contributes to the string tension, and the other 
short monopole loops do not contribute. In the 
other gauges, the monopole density is so dense 
that we can not extract the monopole currents 
which contribute to the confinement. It seems 
that there are many monopole currents which do 
not contribute to the confinement in these gauges. 
We think that these are the reasons why the his- 
tograms in the gauges other than the MA gauge 
do not show so clear change, and that these re- 
sults show that there exist the monopole currents 
which are responsible for the confinement. The 
distributions of these monopole currents change 
at critical temperature in every gauge. 

The calculations were performed on Fujitsu 
VPP500 at the institute of Physical and Chemi- 
cal Research (RIKEN) and National Laboratory 
for High Energy Physics (KEK). 
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